function y = f3(x)

  y = x^4 + x^3 + 3 * x + 1;

endfunction
